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■ Abstract 
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I We use an optimised hopping parameter expansion for the free energy (hnear 

' 5-expansion) to study the phase transitions at finite temperature and finite charge 

density in a global U(l) scalar Higgs sector on the lattice at large lattice couplings. 
We are able to plot out phase diagrams in lattice parameter space and find that the 
standard second-order phase transition with temperature at zero chemical potential 
becomes first order as the chemical potential increases. 
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^ . 1 Introduction 

rS 

' Quantum field theory has had many remarkable successes in the past fifty years using 
the standard methods of perturbative expansions. However, there are many situations 
where such techniques are not appropriate. The obvious example is QCD, which has 
strongly coupled bound states at low momenta. However, phase transitions at non-zero 
temperatures and densities, even in models without infrared slavery like the electroweak 
model, are also poorly explained by perturbative methods. It is this latter example which 
motivates our interest in scalar Higgs fields. The global U(l) model is also a relatively 
simple testing ground for techniques at finite chemical potential and temperature. There 
are many examples of physical situations at finite charge densities. Any phase transitions 
after baryogenesis in the early universe occur at finite baryon charge density. Heavy ion 
collisions also occur in the same regime. 
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Having ruled out a pcrturbativc approach we turn to the various available non- 
perturbative techniques. One's first thought is to use Monte Carlo (MC) techniques 
to tackle the problem. However, MC methods usually fail when considering models at 
finite densities. This is because the Euclidean action, which is used as a statistical weight 
for the system, becomes complex, making a simple statistical integration technique im- 
posiblc. 

To replace the MC approach we need an analytical non-perturbative method. The 
method chosen in this paper is an example of a general family sometimes called lin- 
ear (5-expansions [5]. Examples of these methods have appeared under many names, 
including optimised pertubation theory[4], act ion- variational approach [17], variational 
perturbation theory [21], method of self-similar approximations [22], screened pertur- 
bation theory [20], and the variational cumulant expansion [23]. The method has been 
applied successfully to^ the evaluation of simple integrals [5, 14, 15], solving non-linear 
differential equations [16], quantum mechanics [5, 13, 24, 25, 26], cosmological slow roll 
transitions [7] and to quantum field theory, both in the continuum [26, 21] and on the 
lattice [5, 8, 9, 10, 11, 12, 17, 18, 19, 23, 28, 29]. Since the LDE approach is analytical, 
we do not have to worry about the presence of a complex action. The expectation value 
of all physical observables will turn out to be real. 

The work in this paper with the U{1) or 0{2) model builds on that set out in [4] for 
the case of zero temperature and zero chemical potential. However, here we phrase the 
model in terms of the field and its conjugate ($, $*) rather than working with the real 
components of the field. This change is made because the charge operator is diagonal in 
this representation and so much easier to deal with [32]. 

This paper is also complementary to the work done on the U{1) model using LDE 
methods in the continuum at finite temperature and finite chemical potential [6] . 

2 The Linear Delta expansion 

The general format of the LDE method is to take a given expansion, whether it be a 
perturbative expansion in the continuum or a cumulant expansion on the lattice, and 
to provide an order by order optimisation of this expansion. It is in this process of 
optimisation that non-perturbative information emerges. It is straightforward, in princi- 
ple, to expand beyond leading order, unlike other non-perturbative methods like large-N 
expansions or mean-field approximations. 

The first step in the method is to replace the physical action with an interpolating 
action made up of a linear combination of the physical action and a soluble trial action. 
This trial action is characterised by some set of new variational parameters {v}. The 
particular choice of trial action is the main decision to be made in implementing the 
method. The more general the trial action, and the more terms one includes, the greater 
the number of variational parameters required to characterise it. The interpolating action 
takes the form S Sg — So{v) + 5{S — So{v)). If we set 5 = 1 in the interpolating action 
then Sg — S. The new statistical average of an operator O is 

■^The papers quoted are not necessarily the first in their area but are usually good starting points for 
examining the literature in greater detail. 
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(1) 



We now perform our chosen expansion technique on this new action. In this paper 
we use an expansion for the free energy on the lattice in terms of cumulant averages 
(sometimes known as a hnked cluster expansion). However, we expand in the unphysical 
parameter S instead of the coupling or the hopping parameter. The new expansion, as 
compared to the original, has additional terms at each order which depend upon the 
variational parameters. When we truncate the expansion at a given power of 6 and then 
impose 5—1 there remains a residual dependence on these varational parameters. So 
the LDE up to for the expectation value of some operator O is 



0— »1 



(2) 



where the {}r bracket signifies the truncation of the power series to order S^. It is 
how wc choose to fix the variational parameters, order by order in the expansion, which 
introduces the main fully non-perturbativc effects. There are at least two apparently 
different approaches to this final part of the method listed in the literature. The first is 
to demand that the series converge as fast as possible by minimising the highest order 
term with respect to the lower order terms [17] - the principle of fastest convergence. 
However, we shall take the second main approach, often called the principle of minimal 
sensitivity (PMS)[30]. This is discussed in [5, 8, 9, 10] and used in standard partical 
physics to minimise dependence of perturbative results on the renormalisation scheme 
choice[31]. For the PMS case the parameters, {v}, are set at each order in the expansion 
by demanding that the variation of some physical observable be zero. It is useful to be 
aware of the general statement that a derivative of (2) with respect to some arbitary 
parameter x gives 
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The first term is just an LDE up to order R, the second is just the last term in an LDE 
expansion up to order R and the final term is an LDE up to order R — 1. This means 
that the PMS condition on a variational parameter is 
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In the special case of the free energy, defined by F — —{l/N)\nZ, we have instead 
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In general the PMS condition is only a condition on the final term in the expansion. In 
this way one sees that the PMS and principle of fastest convergence are quite closely 
related. 

The variational parameters are set order by order by (4) or (6) and then substituted 
back into the expression for the physical observable to obtain the LDE estimate for its 
value. In quantum mechanical models at least, the absolute convergence of this LDE 
estimate to the true solution can be proven[13]. 

We have set out the two overall approaches to fixing the varational parameters. How- 
ever, the choice of what physical observable one chooses for the PMS procedure vastly 
increases the possible options. One can for example always apply the PMS techniques to 
an LDE for the particular observable one is trying to measure; this could be called a true 
PMS approach. On the other hand one could choose to extremise an LDE for one partic- 
ular physical observable, say the free energy, and then apply the values of the variational 
parameters gleaned therefrom to another LDE for the particular physical observable one 
is interested in. These two options are discussedin [4] for the 0(2) scalar model and the 
second option is preferred. 

We have set out the two overall approaches to fixing the variational parameters. 
However, there is a further degree of freedom in the choice of which physical observable 
to apply the PMS procedure to. The strict LDE approach would be to apply the PMS 
procedure to the expansion for the observable one is interested in. However, in some 
cases (sec, for example, Ref. 4 for the 0(2) scalar model) there may be no PMS points 
in that expansion. The alternative is to fix the variational parameters by applying the 
PMS procedure to the expansion for one particular, privileged observable, e.g. the free 
energy, and then use those parameters in the expansions for all other observables. 

In this paper we shall take what is effectively the second approach, but instead of 
calculating two expansions, one for the free energy and one for the physical parameter we 
are trying to evaluate, we will evaluate the free energy using an LDE approach with the 
PMS criterion and then calculate all the other required physical observables by numerical 
differentiation with respect to the physical parameters of the free energy. 



3 The U(l) complex scalar model on the lattice 

We work with a Euclidean time formulation of the U{1) field theory. In the continuum the 
global U(l) model at finite density has a conserved charge of the form Q = i($*n — $n*). 
The presence of a finite charge density at finite temperature gives rise to an effective 
action [32] 

-^eff --jy^J [(V^$)*(V''$) + ifxl - ml)^*^ - Ao($*$)' 

-i/Xo(**$ - ***) + Jo* + -^o** (7) 

The temperature appears as a finite boundary condition in the Euclidean time direction, 
where Ntat = 1/T. The integer A^^ is the number of lattice links in the temporal direction 
and at is the temporal lattice spacing. In contrast the chemical potential is present in the 
main body of the Lagrangian. In particular we have a slightly surprising /Xq term which 
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arises from integrating out the conjugate fields when constructing the effective action. 
This term can cause symmetry breaking with positive mp, even for free fields, when the 
transition occurs at /i — mi- In a Euclidean formulation the effective action is complex. 
This complexification caused by the introduction of a chemical potential is generic. In 
addition, the presence of the chemical potential means that the action is not invariant 
under t ^ —t (Osterwalder-Schrader reflection) [1]. 

To regularise the UV infinities of quantum field theory we introduce a lattice formalism 
and replace Euclidean spacetime with a hypercubic 4d lattice. The lattice points are 
denoted by n = (714, n). With a standard redefinition of the physical parameters and 
fields we get 



where Jl — {Ji + %J-i)j\pl. The physical lattice parameters are all dimensionless and 
so measured in units of the lattice cutoff. The derivative terms have become nearest- 
neighbour interactions. We also introduce separate hopping parameters, /t^ and Kt = 1, 
for the spatial and temporal nearest-neighbour interaction. The Kt has a physical value 
of unity but we introduce it as an additional arbitary parameter to allow us to take 
derivatives with respect to it later on. The different hopping parameters allow for a 
different effective lattice spacing in the Euclidean time direction as compared to the 
spatial direction, thus allowing one to indirectly vary the lattice temperature continuously, 
while keeping the number of links in the temporal direction constant, by varing Ks- The 
detail of this process is discussed in section 7.2. To keep calculations relatively easy we 
will keep the Euclidean temporal direction two finks in extent, Nt — 2. This means that 
we will be working with only the first three Matsubara modes and so effectively at high 
temperatures. 

As the physical action stands, when we set = we have automatically, via sym- 
metry arguments, that ($) = 0, so the symmetry is unbroken. This remains the case 
whether one is considering a perturbative expansion in the continuum or a hopping pa- 
rameter expansion on the lattice. 

If the nearest-neighbour terms were to go to zero then the remaining action, containing 
only ultralocal terms, would be exacly soluble. On the lattice we make an expansion 
around this solution in powers of the nearest-neighbour terms, not in powers of the 
coupling. This is often called a hopping parameter expansion. It is the hopping parameter 
expansion for the free energy, often given the special title of a linked cluster expansion 
[35] , which we will optimise in this paper using LDE techniques. 



4 LDE applied to U(l) model on the lattice 

The groundwork for the choice of trial action is set out in [4] . Following the lead of this 
paper we shall include a source term, a quadratic term and a quartic term. Unfortunately, 
because we are dealing with an action in terms of the field and its conjugate rather than 
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its real components, defined by ^ — {(f) + iip) / \/2, we will have to change the notation 
slightly^. 

So = Y^ [{n'' - + Al($"*$")' - - (9) 

n 

where j — {ji + ij^) /\/2. With this choice of 5*0 we obtain an 5'^ which we organise as 
follows: 



Ss = So + 5Si - 5Sint = SuiS) - 6Si 



Sint — ^ ^ 



+Kt{l-llL)^n„nK,+l,n] (10) 

Although the Su is still 6 dependent it is ultralocal, i.e. it contains no interaction 
terms between lattice points. We now perform a diagrammatic expansion of the Si term 
identical to the standard hopping parameter link expansion. As it stands this is not 
strictly a 5 expansion, as there is some residual 5 dependence in the Su action. However, 
this can be dealt with separately later on in the calculation. The expansion is best 
represented graphically, where we define the spatial and temporal links as 



L, := Ks 



n4,a 



n4,n^n4,n+ei 



J]K(l + /iL)$:,,n"^„4+l,n + «;t(l-/XL)$n4A+l,„] := 



S 



int 



(11) 



Thick vertical lines are in the temporal direction, and thin horizontal or diagonal lines 
represent spatial directions. An arrow is added to the temporal link because, due to the 
presence of a non-zero chemical potential, Lj is not symmetric under interchange of initial 
and final lattice points, this is due to the lack of Osterwalder-Schrader symmetry in the 
starting action (7). The arrow can be thought of as representing a net flow of charge in 
the imaginary 'time' direction. 

To extract the thermodynamical information we calculate the free energy per lattice 
site in terms of cumulant averages. To do this we flrst construct statistical averages with 



^To compare the two papers set kj, 



, = ji and jjp = j2 ■ The physical parameters 



are also defined slightly differently: Xejw = 4:Xeim and m^j-^y = 
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respect to the Su action (12); these are given a subscript U. 
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Z = Zu{5){e''')^ = Zu{5)Y, 
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(12) 



Then cumulant averages, denoted by a subscript C, are constructed from the Su statistical 
averages according to 



i=i 



(13) 



We thus obtain the following expression for the free energy 
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(14) 



Because this is in terms of cumulant averages we only have to consider connected dia- 
grams. As an example, the j — 2 term is: 



Et+E 



n,i 
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+2N(d-l){^r)^ + N(d-l){=-)a + N{d-l){. 



(15) 



where d is the number of spacetime dimensions. The multiplicities are calculated by hand 
for each of the distinct diagrams. As one goes up in order the number of diagrams in 
the cumulant expansion increases rapidly. The hopping parameter expansion has been 
evaluated at least up to 14*^^ order for an 0{N) scalar field theory at zero temperature 
and zero chemical potential [33, 34, 38]. There have also been some calculations at finite 
temperature [36, 37]. 

We take the hopping parameter expansion for F up to third order, calculating the 
diagrams by hand. The full list of diagrams and multiplicities used is given in Appendix 
A. Having performed the standard linked cluster expansion to a given order, with the 
additional diagrammatic complications introduced by the presence of the chemical poten- 
tial, we need to convert the expansion to one in terms of powers of 5. To do this we have 
to Taylor expand the cumulant averages to include the residual 5 dependence caused by 
the fact that Su — Su{S). This will introduce additional optimisable terms at each order 
in the expansion. This method of caknilating the LDE is much simpler than including all 
the S dependent terms in Sint at the start, which requires many additional diagrams. 

Not all the diagrams in the hopping parameter expansion 'feel' the temperature, as 
they are not periodic across the Euclidean time direction. However, when we apply the 
variational part of the LDE approach all the diagrams effects are combined and so the 
temperature is felt indirectly by the whole of the expansion. 
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The presence of a j*$ term in the trial action, where the j value is fixed variationally, 
allows for the possibility that ($) 7^ even when Jl = 0. This means that we can 
examine the symmetry breaking of the model using the optimised LDE version, unlike 
the standard free energy expansion, which is fixed on one side of the phase transition by 
the symmetry of the physical action. 



5 Example diagram calculation 

The statistical and cumulant averages result in an expression made up of known param- 
eters and a set of integrals involving the Lu Lagrangian. These integrals in turn, when 
one performs the Taylor expansion in 5, become integrals with respect to Lq because 
Lu — Lq + 5Li. At this point we should define the general set of integrals 



JmniS) := J d(j)d(P*{(Pr{(p*re-''^^'^ Kmn{S) := J #2#l(0l)"X02)"e 
Bmn := j #2#l(0l)"*(02)"e-^° 



-Lu{5) 



Joo{S) Koo{d) Boo 



In terms of these definitions. 



lnZ[/ = 7VlnJoo (17) 

All the diagrams we need to calculate are expressible in terms of these integrals. 
Consider the very simple example of the cumulant average of the single spatial link. The 
cumulant as defined in (13) is expandable in terms of statistical averages with respect to 
the Su action (12). We get 



( )c = ( ) 
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n4,nC4,n+ei>t7 = l^s [{<t>*a4>b)u + {4>a4'l)u\ 



/(a 



/(n„.#mrf0^)exp{-En^&(<^)} 



(18) 



As Lu is ultralocal all the integrals separate from each other. It was to ensure this 
fact that we imposed ultralocality on our choice of trial action, 5*0. The integrals over 
variables other than 0^ and 0b cancel between the numerator and denominator. Therefore 
in both terms in (18) we will just be left with the product of two integrals. Using the 
identities defined in (16) we get 
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(19) 



All the cumulants obtained from the expansion in (14) are similarly expressible in 
terms of some combination of our integrals. The final step in evaluating the expansion of 
F to a given explicit order in 6 is to Taylor expand the J integrals. We first express the 
j integrals as K integrals. Then we Taylor expand the later using the following recursion 
relation for the derivative: 
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Noting that 



(5=0 



Bmn we have all the required information to evaluate the 

Taylor expansion to any order. For example, the single spatial link can now be calculated 
using 
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(21) 



(22) 



All the other cumulant averages, initially evaluated in terms of J integrals, and the 
In Joo term, can be similarly expressed as a power series in 5 in terms of the B integrals. 
As one increases the order of the 5-expansion the number of diagrams soon becomes too 
large to deal with by hand. At order 3, for example, there are roughly 20 diagrams (see 
Appendix A). We therefore use an algebraic manipulation program which can handle 
the repetitious task of calculating all the Taylor expansions and rearranging the overall 
expansion in explicit powers of 5. 



6 Minimisation 

As the expression for F given in (14) stands, each order in the expansion has residual 
dependence on the variational parameters {fi^, ji, ^2}, as well as the expected dependence 
on the physical parameters {m|, k^, kj, /x^,, Al, Ji, J2}. We will fix the variational param- 
eters using the PMS condition set out in (6). But first we note that we can immediately 
set j2 = 0, which can be justified using symmetry considerations [4]. Having made this 
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choice we are reduced to a two-dimensional variational parameter space. Minimisation 
of the free energy in this space is best represented as a two-dimensional contour plot. 
The phase transitions in the physical parameter space are actually mirrored in the vari- 
ational parameter space in terms of the number of minima present. This behaviour has 
already been noted in more limited terms in [4]. As a result of this there are a few generic 
pictures for the minimisation depending on the values of the physical parameters. The 
minimisation curves are plotted in Figure 1. 

We include the lattice temperature with these plots. See section 7.2 for the derivation 
of these temperatures. Going from Fig. 1(a) to 1(b) we increase Kg and so decrease 
the lattice temperature. In so doing we move from one minimum to two minima in the 
variational parameter space and at the same time undergo a second order phase transition 
in the physical parameter space. Similarly, going from Fig. 1(a) to 1(c) we increase hl 
(also indirectly increase T) and in so doing move from one minimum to three minima 
in the variational parameter space. At around the same time in the physical parameter 
space we undergo a first-order phase transition. The order of the 'phase transition' in 
variational parameter space matches the order of the true phase transition in the physical 
parameter space. 

In the unbroken physical phase the minima will be found at the point: 

q2 ^ ji = (23) 

This is discussed further for the zero temperature and zero chemical potential case in [4] . 
As noted in that paper, because the minimisation space for the free energy reflects the 
physical parameter phase space, the value of the variational parameters at the minimum 
is often sufficient information to find the phase of the system, without a full calculation 
of ($). However, in the case of a first-order phase transition more care is needed, as 
there will be both 'broken' and 'unbroken' minima, as in Fig. 1(c). One chooses between 
them by taking the overall minimum. As an aside, note that minimising the free energy 
is equivalent to maximising the entropy, which gives us a further physical justification for 
our particular PMS approach. 

By tracking the value of F at all the minima, and taking special care when we have 
a first-order transition, we can build up a set of F values for any range of physical 
parameters. From the free energy data we can then evaluate many other correlators by 
taking numerical derivatives with respect to the physical parameters. 

Of course the method is always limited by the fact that it requires the presence of 
a minimum. When the minima in F disappear the method breaks down. However, the 
method is sustained far enough on either side of the phase transition for a large enough 
range of physical parameters to allow it to be used to plot out reasonable phase diagrams. 

Having minimised the free energy we can now take numerical derivatives of F to 
evaluate particular physical observables. This will allow us to build up phase portraits 
in the {T, n^} plane for any choice of the other physical parameters. 
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One minimum at ~ 
Two saddle points at fi^ ~ —Ax 

Ks - 0.6, /YL - n,T - 2.5 



Two minima at ~ m| 
Two saddle points at fi^ ~ — 

K<j = 1, /Ui = fl, 7" = fl-77 




One minimum at fi^ ~ 
Two minima and two saddle points at 

K.S — 0.6, = ?i^T = 5.93 



One saddle point at Q? ~ 
-Ax Two minima and two saddle points at fi^ ~ —Ax 

= 1, = 'i-,T = 4.33 



Figure 1: Various minimisation plots. In all graphs m\ — — 40, = 1, Al = 100, Ji = 
J2 = 0. 
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7 Numerical Results 



7.1 Phase transition using Kg 

As a first example we shall look at evaluating as a function of the hopping parameter, 
Ks, to look for a phase transition. To do this we evaluate the estimated F value, by 
tracking the minima across our 2D variational parameter space, at two nearby Ji values 
around zero to allow us to evaluate a numerical derivative with respect to Ji. This is 
related to the value using (5) to give 




Therefore by taking numerical derivatives we can look for symmetry breaking, where 
7^ 0. For simplicity we consider the /i^ = curves. In this case the variational 
parameter space tracks between Fig. 1(a) and Fig. 1(b). 




Figure 2: Phase transition with Kg 



We have a second order phase transition to a broken symmetry sector as we increase 
the Kg value, and so decrease the lattice temperature. This transition happens at lower 
and lower temperatures as one increases the lattice mass. 

7.2 Evaluating the lattice temperature 

To construct a true phase diagram we need to translate from the hopping parameter Kg to 
the lattice temperature T. Therefore we need to calculate the lattice temperature caused 
by particular choices of the Kg hopping parameter. 
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To do this wc note that the presence of different hopping parameters in the tempo- 
ral and spatial directions means that correlation functions will also be anisotropic with 
respect to the temporal and spatial axes of the lattice. The correlation length in lattice 
units in the time direction, ^t, will be different from the space-like correlation length, ^g- 
However, we expect the correlation lengths in physical units to be the same, even though 
the lattice spacings at and will be different. This means that a^^s = at^t, which implies 
that 

| = ^ = 2Ta, = ^ (25) 

where is the spatial cutoff. By varying the spatial hopping parameter. Kg, we will vary 
the relative size of the two correlation lengths and so vary the effective lattice anisotropy. 
This in turn allows us to vary the temperature [1]. 

The overall scale of the lattice temperature is set by A^t = 2, which hmits us to 
high temperatures, but within this range of high values we can vary the temperature 
continuously using the hopping parameter Kg- 

Now (5) gives 



dF 



= L / ^ \ ^ -(^d _ 1) ^^n4,n*^n4,n+ei _^ ^n4,n^n4,n+ei*^ ^ _ j^-j^^ 




OKs N 

dKt ^ \ n I 

(d-lldF /dF 
As 2 dKt / dKs ^ ' 

Therefore purely by taking numerical derivatives with respect to the hopping param- 
eters we can construct the lattice temperature for a given value of Kg. As an example we 
plot ln(T/A5) against In^Kg) for a particular range of physical parameters in Figure 3. 

In the unbroken regime of physical parameters there is a simple inverse relationship 
KgT = k. This is true for both zero and non-zero chemical potential. The only effect of 
the chemical potential is to change the value of the constant k. In the unbroken regime 
we find an almost perfect fit for the empirical curve 

KsT=^{fil + 3) (27) 

However, as we pass through what is a second order phase transition, the relationship 
becomes more complex. This break-off point occurs at lower and lower Kg as one decreases 
the lattice mass. 

For any other set of physical parameters we can similarly evaluate the effective tem- 
perature given by a particular choice of Kg. 

7.3 Phase transition using T 

We can now combine the information in sections 7.1 and 7.2 to allow us to evaluate ($1) 
as a function of T and look for a phase transition. This leads to the curves shown in 
Fig. 4. 
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ln(r/A,) 




Kt = 1,Al = 100, Ji = 0,^2 = 



Figure 3: Effective T 




Figure 4: Phase transition with T 
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These curves follow directly from Fig. 2 and, as expected, as the temperature is 
lowered the symmetry breaks at some point, starting first with the lowest mass. The 
phase transition is second order in form. As the temperature decreases towards zero 
there is an unexpected divergence of the value of However, this is precisely in the 

regime where the approximation we are using is expected to break down, so no physical 
interpretation should be made of this behaviour. 

7.4 Phase transition using m\ 

For completeness we show the standard symmetry breaking which occurs as one lowers 
the physical mass. This second-order phase transition is discussed in much more detail 
in [4]. For particular choices of physical parameters we have the phase transition curves 
shown in Figure 5. 




Kt = l,llL = 0,XL = 100, Ji = 0, J2 = 

Figure 5: Phase transition with m| 

As expected, lowering the lattice mass breaks the symmetry, starting first with the 
highest Ks curve. The phase transition is second order. The curves are performed at fixed 
hopping parameter values and not at fixed lattice temperature. In the unbroken phase, 
as we can see in Fig. 3, the effective lattice temperature is independent of the lattice 
mass, so fixing Kg is equivalent to fixing T. However, when the symmetry is broken the 
temperature becomes mass dependent. Although this make things more complex it does 
not stop us from measuring the lattice temperature at the phase transition. Using the 
empirical formula (27) we can calculate the lattice temperatures. Ks = 0.6,0.7,0.8,0.9 
and 1.0 correspond to T = 2.499, 2.143, 1.875, 1.667 and 1.5 respectively. 
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7.5 Phase transition using fi^ 



Having seen how the symmetry is restored as one increases the temperature and lattice 
mass, we should finally look at evaluating as a function of fii- The phase transition 
with involves travelling between variational parameter spaces which look like Fig. 1(a) 
and Fig. 1(c). This means that we will have first-order phase transitions. Examining the 
free energy across the transition leads to the curves shown in Figure 6. 




Figure 6: Phase transition in F with fi^ 



In the m| = —15 • ■ ■ — 25 plots we can see that there is a kink in the free energy 
due to the first order phase transition from an unbroken (u) to a broken (b) global U{1) 
symmetry. The transition point moves to lower /i values as the mass is decreases. These 
curves have variational parameter spaces which move between Fig. 1(a) and Fig. 1(c). 
The kinks arise as one jumps from the unbroken to the broken minima in the variational 
parameter space. 

At m| < —30, and lower, the symmetry is broken for all fi values. The variational 
parameter space moves between Fig. 1(b) and Fig. 1(d). The symmetry is already broken 
ai = because of the second order phase transition which occurs at low enough 
temperatures, as discussed in section 7.4. The special case of the m| = —30 curve is 
discussed in more detail later. 

Using the free energy data, and again noting (24), we can construct the value. 
Where there is a kink in the F value we expect a first order phase transition in the field 
value. This is precisely the behaviour seen for the curves in Figure 7. 

In the m| = — 15 ■ ■ ■ — 25 plots we see there is a first order phase transition from 
an unbroken to a broken global U{1) symmetry. Crossing through this phase transition 
we move from Fig. 1(a) to Fig. 1(c) in variational space. Raising the chemical potential 
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Figure 7: Phase transition in with fj,L 



breaks the symmetry, starting first with the lowest curve. The transition point moves 
to higher fi^ values as one raises the lattice mass. 

The m| = —30 case is more unusual because, although the symmetry is already 
broken at /i^ = 0, we find evidence for a further phase transition at /i^ = 2.1. It turns 
out that this unexpected behaviour occurs at the place in phase space where the first order 
phase transition meets the second order phase transition. We shall postpone discussion 
of this point, as it will become much clearer on consideration of the full phase diagram, 
see section 7.6. 

For m| < —30 we see that the symmetry is already completely broken. There is a 
smooth crossover as one increases the chemical potential value from zero and in doing 
this one moves from Fig. 1(b) to Fig. 1(d) in variational parameter space. 

Just as in section 7.4 these curves are derived at fixed hopping parameter values 
and not at fixed lattice temperature. However, unlike the lattice mass, the temperature 
is chemical potential dependent, whether or not one is in the broken regime. Each 
point on the curves is at a different lattice temperature. Again this make things more 
complex but we can still find the lattice temperature at the phase transition. We find 



that that {ml = -15, /i^ = 4} gives T = 9.5, {ml = -20, /i^ = 3.5} gives T = 7.625, 
{ml = -25, /iL = 2.9} gives T = 5.705, and finally {ml = -30, /il = 2.1} gives 



7.6 Phase diagrams 
7.6.1 {T,hl} Phase diagram 

Combining all the information from the previous sections we can build up phase diagrams 
for the theory in {T, fi^} space. These are plotted for various value of the lattice mass in 



T = 3.705. 
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Figure 8. 
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Figure 8: {T, /i^} Phase diagram 

Close to the iil = axis we have a second-order phase transition between the unbroken 
phase at high temperatures and the broken phase at low temperatures. As one crosses 
through this transition the minima plots change from Fig. 1(a) to Fig. 1(c). However, as 
one follows this transition out to higher densities it becomes first order in nature. For 
the first order phase transition the minima plots change from Fig. 1(a) to Fig. 1(b). 

For a small region at the changeover between the first-order and second-order phase 
transition lines both the minima at — A^, and the minima at m| are present in variational 
parameter space. Travelling from the second-order regime into the first-order part the 
minima at m| spawn two further minima which rapidly move up towards — Al- This 
is captured in the variational parameter space picture seen in Figure 9. 

As fiL and T are increased further one enters the first-order phase transition regime 
and the two minima at m| become one minimum. 

It is this overlap which explains the unusual phase structure seen for the m| = —30 
curve in Fig. 7. We can that the unexpected additional first-order transition for this 
curve sits at fii = 2.1,T = 3.7, which is right in the changeover regime for the m| = —30 
case in Fig. 8. For fiL < 2.1 the field value is being evaluated at one of a pair of minima 
at m|, for the fii > 2.1 the field value is being evaluated at a pair of minima at — A^. 
The jump between the two causes the first-order phase transition. 

In physical space we have the three distinct types of pictures on the broken side of 
the transition, these are shown in Fig. 10. 

As one moves through the changeover area, starting from the second-order area of the 
phase diagram, the first-order style metastable states appear out beyond the second-order 
style minima. Then as one moves into the first-order area the two second-order minima 
merge to give the classic first-order field distribution. 
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Figure 9: Changeover region 



7.6.2 {T,p} Phase diagram 

One can also calculate the charge density in lattice units which, noting that Tie// 
7i + /iQ and = at^, is 



P ■-- 



iQ) 

Nsal 



2T 



dF 



{2i 



where, because we are in lattice units, = = 1. Using (28) one can construct the 
phase diagram of temperature against charge density. This is plotted in Figure 11. 

The lower ends of the transition lines move closer to the origin in {T, p} space as one 
increases the mass. In the close-up of the second-order phase transition lines, near to the 
/X = axis, we can see an interesting complex curve form. As the transition becomes 






2nd order 



changeover 



1st order 



Figure 10: Qualitative F against ($1) curves along the phase transition line 
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Figure 11: {T,p} Phase diagram 
first-order for higher chemical potentials the transition lines become straighter. 

8 Conclusions 

An LDE optimisation of the standard hopping parameter expansion has allowed access 
to some of the truly non-perturbative physics of the scalar global U{1) model. Close to 
the /ii = axis we find the expected second-order phase transition but as we track out to 
higher chemical potentials we see that this transition becomes first order. At our level of 
approximation at least there appears to be evidence for a more complex phase structure 
than expected. 

Although the U{1) model is in itself intrinsically interesting as a simple model for 
Higgs physics, the work in this paper can also be seen as a test of LDE techniques as 
applied to phase transition at finite density in general. In this context we can see that 
the LDE is a sucessful approach. The presence of the additional variational parameters 
allows one to examine a theory as one tracks through the phase transition with any of 
the physical parameters. By the choice of variational parameter one can actually choose 
which particular equivalent broken vacuum state the model is in, as demonstrated in 
[4]. Also the form of the free energy contours in variational parameter space signals the 
phase-space behaviour of the physical model itself. This allows one to extract the physical 
phase of the system without explicitly calculating the particular order parameter such as 

To extend our approach to plotting out a phase diagram at finite density to a more 
complex model, in particular a gauge theory, we need to find an equivalent, lattice reg- 
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ulariscd, expansion to optimise. The hopping parameter expansion is only available for 
scalar theories: the derivative terms cannot be broken up in the same way for a gauge 
field theory. However, there is an alternative lattice expansion in the strong coupling 
expansion, i.e. an expansion in powers of the trace round lattice plaquettes. An example 
of an LDE optimisation of a plaquette expansion for gauge theories on the lattice at zero 
temperature and chemical potential is given in [8]. Using a similar scheme to that used 
in this paper for scalars one could extend considerations of gauge theories on the lattice 
to finite temperature and densities. 
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A Diagrams 

Below we list all the diagrams which contribute to the cumulant expansion for the free 
energy of a U{1) scalar Higgs model with finite temperature and chemical potential up 
and including order 3. With a lattice of two temporal links in extent we have periodicity 
across some of the diagrams. Points which are identified as identical are signified with an 
open circle in the diagrams. In the multiplicities that accompany the diagrams we have 
omitted the overall N factor and used the shorthand 

dj = 2{d-l)-j, (29) 

where d is the spacetime dimension of the lattice and j & Z. 
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m2,3 = 2do 



m2,5 = 1 
"^3,1 = Ido 
m3,3 = ^dodl 

"^3,7 = Qdodi 
7713,9 = 6dg 



^3,11 = 6dodi 

^3,13 = 3do 



mi,2 = 1 
m2,2 = dodi 

777,2,4 = 2do 



7712,6 = 2 
"^3,2 = 3dodi 
77l3,4 = (io(ii(i2 

"73,6 = 3(io 
7773,8 = 3cioC?l 
?773,10 = 3doC^l 



r773,i2 = Mo 
7773,14 = 3do 



n 



?773,17 = 3do 



?773,19 = 3 



?773,16 = 3do 



^3,18 = 1 
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